SMOOTH APPROXIMATION OF BI-LIPSCHITZ 
ORIENTATION-PRESERVING HOMEOMORPHISMS 



SARA DANERI AND ALDO PRATELLI 

Abstract. We show that a planar bi-Lipschitz orientation-preserving homeomorphism can 
be approximated in the W 1,v norm, together with its inverse, with an orientation-preserving 
homeomorphism which is piecewise affine or smooth. 



1. Introduction 

In this paper we deal with approximations of bi-Lipschitz orientation-preserving homeomor- 
phisms u : f2 C M 2 — > ACM 2 , where $7 and A = u(Cl) are two open bounded subsets of 
M 2 . In particular, we show that both u and its inverse can be approximated in the ' p -norm 
{p € [1, +oo)) by piecewise affine or smooth orientation-preserving homeomorphisms. Our main 
theorem is the following. 

Theorem 1.1. Let Q C K 2 be any bounded open set, and let u : Q — > A be a bi-Lipschitz 
orientation-preserving homeomorphism. Then, for any e > and any 1 < p < oo, there exists a 
bi-Lipschitz orientation-preserving homeomorphism v : Q — > A such that u = v on dfl, 

lh ~ v \\l°°(Q) + Wu' 1 ~ IIl°°(A) + \\Du - Dv\\ LP{n) + \\Du~ 1 - £>u -1 [|.lp(a) < e, (1.1) 

and v is either countably piecewise affine or smooth. More precisely, there exist two geometric 
constants C\ and C2 such that, if u is L bi-Lipschitz, then the countably piecewise affine approx- 
imation can be chosen to be C±L 4 bi-Lipschitz, while the smooth approximation can be chosen 
to be C2L 28 / 3 bi-Lipschitz. 



Thanks to a result by Mora-Corral and the second author [28J (see Theorem 2.1 below), the 



problem of finding smooth approximations can be actually reduced to find countably piecewise 



affine ones -i.e. affine on the elements of a locally finite triangulation of f2, see Definition 3.2 

The fact that v might not be (finitely) piecewise affine but countably piecewise affine is due 
to the fact that we require u = v on d£l, so it is clearly impossible to find a (finitely) piecewise 
affine approximation v unless the domain is a polygon and u is piecewise affine on the boundary. 
In fact, we also prove the following result. 

Theorem 1.2. If under the assumptions of Theorem \l.l\ one has also that £1 is polygonal and u 
is piecewise affine on dQ, then there exists a (finitely) piecewise affine approximation v : — > A 
as in Theorem \l.l\ which is C\C'{£l)L* bi-Lipschitz. 



About the dependence of C '(O) in Theorem 1.2 on the domain f2, see Remark 6.1 
The first naive idea coming to one's mind in order to construct a piecewise affine approxima- 
tion of u could be the following: first, to select an arbitrary locally affine triangulation of with 
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triangles of sufficiently small diameter; then, to define v as the function which, on every triangle, 
is the affine interpolation of the values of u on its vertices. Unfortunately, if on one hand the 
functions defined in this way provide an approximation of u in L°°, on the other hand they may 
fail to be homeomorphisms. The problem is due to the fact that, taking arbitrary nondegenerate 
triangles in Q -no matter how small- then the affine interpolation of u on the vertices of the 
triangles can be orientation-preserving on some triangles and orientation-reversing on the others 
(see Figure [I]). This prevents the affine interpolation to be injective since an homeomorphism 
on a connected domain in IR 2 must be either orientation-preserving on every subdomain, or 
orientation-reversing on every subdomain. An explicit example of a function with such a bad 
behaviour can be found in [3"T] . 




FIGURE 1. The square ABCD is divided in the triangles T and T'. The affine interpo- 
lation v of u on ABCD is not injective, since v(T) C v(T') (v(T) and v(T') are shaded). 
Moreover, u is orientation-preserving in the square while v is orientation-reversing on T . 

The general problem of finding suitable approximations of homeomorphisms u : M rf I) ^ — > 
u(fi) C R d with piecewise affine homeomorphisms has a long history. As far as we know, in 
the simplest non-trivial setting (i.e. d = 2, approximations in the L°°-norm) the problem was 
solved by Rado [29]. Due to its fundamental importance in geometric topology, the problem of 
finding piecewise affine homeomorphic approximations in the L°°-norm and dimensions d > 2 
was deeply investigated in the 50s and 60s. In particular, it was solved by Moise [25] and Bing [3] 
in the case d = 3 (see also the survey book j2H] ) , while for contractible spaces of dimension d > 5 
the result follows from theorems of Connell [8] , Bing [6] , Kirby |22j and Kirby, Siebenmann and 
Wall [23] (for a proof see, e.g., Rushing [30J or Luukkainen [24J). Finally, twenty years later, 
while studying the class of quasi-conformal varietes, Donaldson and Sullivan [12] proved that 
the result is false in dimension 4. 

Let us now consider the case of Sobolev homeomorphisms u £ W 1,p for some p G [l,+oo]. 
As pointed out by Ball (see [21 [3] , see also Evans [13] ) , the problem of proving the existence of 
piecewise affine approximations of Sobolev homeomorphisms arises naturally when one wants to 
approximate with finite elements the solutions of minimization problems in nonlinear elasticity 
(e.g. the minima of neohookean functionals, see also pQ, [1], [9], [32]). In that context, the 
function u represents the physical deformation of a material with no interpenetration of matter 
(in particular, d = 2 as in the present paper, or d = 3). 
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The results available in the literature provide, under increasingly weaker hypotheses on the 
derivatives of u, piecewise afhne or smooth approximations of u and its derivatives. The first 
results were obtained by Mora-Corral [27J (for planar bi-Sobolev mappings that are smooth 
outside a finite set) and by Bellido and Mora-Corral [7J, in which they prove that if u E C 0,a for 
some a E (0, 1], then one can find piecewise affine approximations v in C ' 13 , where f3 E (0, a) 
depends only on a. 

Recently, Iwaniec, Kovalev and Onninen [21] almost completely solved the approximation 
problem of planar Sobolev homeomorphisms, proving that whenever u belongs to W l ' p for some 
1 < p < +oo, then it can be approximated by smooth diffeomorphisms v in the W 1,p -noim 
(improving the previous result for homeomorphisms in W 1,2 found in [20j). 

In all the above-mentioned results, approximations are intended in the W 1,p sense, that is, 



in place of (1.1) one obtains estimates of the form 

||« - v\\lp(Q) + \\Du - Dv\\ LV (si) < e, (1-2) 

whitout any information on the inverses u~ l and v~ 1 . Unfortunately, this is still not enough for 
the applications: in fact, the functionals of nonlinear elasticity usually depend on functions of 
the Jacobian of u which explode when det(Du) — > (see for instance [31 pag. 3]). The physical 
meaning of choosing such functionals is that too high compressions or strecthings require high 
energy. As a consequence of this, two invertible Sobolev functions u and v which are close in 



the sense of (1.2) may have very different energies. Instead, if u and v are close in the sense 
of (1.1), then their energies are also close. 

These considerations suggest to work in the space of bi-Sobolev homeomorphisms, that is, 
homeomorphisms u such that both u and u _1 belong to W 1 '?. This was also already suggested 
in the paper by Iwaniec, Kovalev and Onninen ( |21| Question 4.2]). We only mention here that 
the study of bi-Sobolev homeomorphisms is very active, also in connection with maps of finite 
distorsion (see for instance [M El ESI E3 13 HSl IS]). 

The present paper is the first one to take care also of the distance of the inverse maps, 
leading to a partial result towards the solution of the general problem (hence, we also give a 
partial positive answer to Question 4.2 of |21j): in fact, we are able to deal with homeomorphisms 
which are bi-Sobolev for p = +oo. The techniques adopted in [7J and [21] are completely different 
with respect to the ones which will be used throughout this paper. While the proof in [7J is 
based on a refinement of the supremum norm approximation of Moise |25] (which, as pointed out 
by the authors themselves, cannot be extended to deal with the Sobolev case) and the approach 
of |21j makes use of the identification M? ~ C and involves coordinate- wise p-harmonic functions, 
our proof is constructive and based on an explicit subdivision of the domain of u depending on 
the Lebesgue points of Du. 

We conclude giving a bound on the values of the costants C\, C2 and C3 appearing in 



Theorems 1 1 . 1 1 [L~2~| and 2.2 (while the constant C'{Q) depends on the set Q, see Remark 6.1), 



Ci = 72 4 C 3 , C 2 = 70Cl /3 , C 3 = 636000 . 
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2. Scheme of the proof and plan of the paper 

Our proof is constructive, thus long, but it relies essentially on three known facts: the 
Lebesgue differentiation Theorem for L 1 -maps in M d , the Jordan curve theorem and a planar 
bi-Lipschitz extension theorem for homeomorphic images of squares proved in [llj . As already 
mentioned in the introduction, all our effort will be to get a piecewise affine approximation of 
u, since then the smooth extension readily follows by the following recent result from [28J. 

Theorem 2.1. Let v : SI — > IR 2 be a (countably) piecewise affine homeomorphism, bi-Lipschitz 
with constant L. Then there exists a smooth diffeomorphism v : O — > u(Q) such that v = v on 
dVL, v is bi-Lipschitz with constant at most 70L 7 / 3 , and 

\\v - w|U«(n) + \\Dv - Dv\\up(q) + - v -1 ||i,«>( v (n)) + \\Dv~ 1 - Dv^Wlp^^ < e . 

A rough idea and scheme of our construction is as follows. 
Approximation of u on Lebesgue squares. The first idea is to use the fact that, in a 
sufficiently small neighborhood of each Lebesgue point z for the differential Du, the map u 
is arbitrarily close, both in W ' p and in L°° , to an affine L bi-Lipschitz map (given by its 
linearization around the point z). The W ,p estimate is simply a restatement of the definition 



of Lebesgue point of Du, while the L°° estimate is proven in Lemma 4.2 Indeed we prove 
that, given a square DC!1 (e.g. a neighborhood of z), the more Du is close in L P (D) to an L 
bi-Lipschitz matrix M (given e.g. by Du(z)), the more u is close in L°°(T)) to an L bi-Lipschitz 
affine map um with Dum = M. Moreover, since u is bi-Lipschitz, we have that also the inverse 
of u is close both in W 1,p (u(V)) and in L°°(u(V)) to the inverse of um- 

The main implication of these estimates towards the construction of a piecewise affine bi- 
Lipschitz map approximating u is the following. Let us take a square PCfias above and let 
us consider the piecewise affine function v which coincides with u on the vertices of T> and is 
affine on each of the two triangles obtained dividing T> with a diagonal. If \\Du — M\\jj,(p\ is 
sufficiently small, then the L°° estimate implies that u(dT>) is uniformly relatively close to the 
parallelogram of side lengths at least side(P)/L given by um{9V). Hence, since v = u on the 
vertices of the square and is affine on each side of dT>, the same uniform estimate holds also 
for v. In particular, the map v is orientation preserving, injective, and approximates u and its 
inverse as desired. 

Finally, thanks to the fact that the Lebesgue points of Du have full measure in O, we fix 
two orthonormal vectors ei, e2 £ K 2 and, Ve > 0, we find a set VL £ CC SI with j£f ($7 \ Q e ) < e 
which is made by a uniform "tiling" of squares with sides parallel to ei, e2 with the following 
property. On each square T> of the tiling, Du is sufficiently close to an L bi-Lipschitz matrix M 
(in particular, M will be equal to Du{z) for some Lebesgue point z € D). Then, by the previous 
remarks, one can show that the piecewise affine function v obtained interpolating between the 



values of u on the vertices of the squares is injective and satisfies (1.1) on $7 e . Moreover, v is 
L + e bi-Lipschitz. The squares of the tiling covering Vt £ will be called Lebesgue squares, and the 
set J\ right polygon, due to its shape -see Figure [2j 
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Thus, the first idea of the proof leads to define a piecewise affine approximation of u on a set 
whose Lebesgue measure is as close as we want to Jzf (f2). In order to complete the construction, 
we have to define v in the interior of the set Q, \ Q £ . 

Countably piecewise affine bi-Lipschitz extension. The second idea of our proof is to 
reduce to the following model case: $7 \ S\ is a square of Lebesgue measure at most e and 
n la(n\n ) * s a piecewise affine function. In particular, by the previous construction, v = u on 
d£l £ . In this case, an approximating function v is provided by the following bi-Lipschitz extension 
theorem, proved in 



Theorem 2.2 ([IT]). There exists a geometric constant C3 such that every L bi-Lipschitz piece- 
wise affine map u : &D(0, 1) — > M 2 defined on the boundary of the unit square admits a C3L 4 
bi-Lipschitz piecewise affine extension v : V(0, 1) — > T, where T is the bounded closed set such 
that dT = u(dV(0,l)). 

In particular, it is shown in that one can take C3 = 636000. We recall that an analogous 
result had already been proved by Tukia in [33]. However, explicit estimates of the Lipschitz 
constant of the extension v were not provided. 

It is then easy to verify that, provided e is chosen sufficiently small at the beginning, such an 
extension of u^ g ^ n } together with the already defined piecewise affine interpolation of u on the 
Lebesgue squares, satisfies the assumptions of Theorem Indeed, by definition, v is injective 



on the whole Q. Moreover, we know by the previous construction that it satisfies (1.1) on Q £ . 
And on the other hand, on Q \ £l £ , \Du\ and \Dv\ are bounded by the two Lipschitz constants 
L and CL 4 (together with their inverses) on a set of small area and then the W 1,p estimates 



in (1.1) follow. Finally, since Q\f2 e and u(0\J7 E ) have small Lebesgue measure, v and v^ 1 are 
also close to u and in L°° . 

In order to reduce to this model case, we perform the following steps: 



6 



SARA DANERI AND ALDO PRATELLI 




1. We cover Q\Q e with a countable (locally finite in $7) "tiling" of small squares whose 
sides are parallel to ei and e2 (see Figure [3]). 

2. On the 1-dimensional grid Q given by the boundaries of the squares of the tiling we define 
the piecewise afhne approximation v in such a way that v(Q) C A and v is 72L bi-Lipschitz. 



3. We "fill" the squares of the tiling extending v\ by means of Theorem 2.2 thus getting 
a globally C*3(72L) 4 bi-Lipschitz function on £2 \ fl £ . 

The fact that the Lipschitz constant of v on £l e depends only on the Lipschitz constant of 
u will tell us that, as in the model case, the W 1,p and L°° norms of u — v and u _1 — v~ l can be 
made as small as we want -provided we choose e sufficiently small at the beginning. Thus we 
end the proof of Theorem 

Let us also give a very rough idea of how the proofs of Steps 1, 2 and 3 works. While 
Step 1 is a simple geometric construction, Step 2 essentially consists in approximating u on the 



grid Q with a piecewise affine function. This will be possible thanks to Lemma 5.5, which tells 



that it is possible to approximate u on the segments, and to Lemma 5.6, which takes care of the 
"crosses". Finally, Proposition 5T in Section [5] concludes the argument of Steps 2 and 3. The 
essential idea there is that, since on the "non-Lebesgue squares" the behaviour of u is wilder, 
one cannot simply take v equal to the affine interpolation of u on the vertices. Indeed, as already 
pointed out in the introduction, this could easily give a non-injective function. However, since 
the total area of the non-Lebesgue squares is small, any bi-Lipschitz function which preserves 
the boundary values approximates u as in (1.1). This is why we use the extension of v\q given 
by Theorem [272] in Step 3. 



Plan of the paper. Section [3] contains the main notation and some preliminary definitions. 

Section [4] is devoted to the construction of v on a large right polygon made of "Lebesgue 
squares" . 
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In Section [5] we complete the proof of Theorem 1.1 using the reduction argument outlined 
above and the bi-Lipschitz extension Theorem |2.2| 



Finally, in Section [6j we adapt the proof of Theorem |1.1| to provide a (finitely) piecewise 
affine approximation of u whenever Q is polygonal and u\ gn is piecewise affine, thus proving 
Theorem 11.21 

3. Preliminary Notation and Definitions 

In this section we give some preliminary definitions and fix some useful notation which will 
be used throughout the paper. 

First we recall the definition of orientation-preserving (resp. reversing) homeomorphism. 

Definition 3.1 (Orientation-preserving (reversing) homeomorphism). We say that an homeo- 
morphism u : Q — > u(Q) C IR 2 is orientation-preserving (reversing) if whenever a simple closed 
curve [0, 1] 3 t t-t ^y(t) G Q, is parameterized clockwise, then [0, 1] 3 t i-> u(7(i)) G u(Q) is 
parameterized clockwise (resp. anti-clockwise). 

It is well known that if $7 is connected, then any homeomorphism u : Cl — > u(O) C M 2 is 
either orientation-preserving or orientation-reversing. 

Next, we define the class of functions in which we look for approximations of u. To this aim 
we recall the definitions of (finite ) triangulation of a polygon and of locally finite triangulation 
of an open set fi C IR 2 . 

Definition 3.2 ((Finite) triangulation). A (finite) triangulation of a polygon Q,' CM? is a finite 
collection of closed triangles {Ti}?L 1 whose union is equal to closfi' and, for all i 7^ j, 

T{ n Tj is either empty, or a common vertex, or a common side ofTi and Tj . (3.1) 

Definition 3.3 (Locally finite triangulation). Let Q, C IR 2 be a bounded open set. A locally 
finite triangulation of Q is a locally finite collection of closed triangles {Tj}j g N such that ft C 



UieN-^ — dos^ an d (3.1) holds. 



We notice that, unless fl is a polygon, the number of elements of a triangulation cannot be 
finite. 

Definition 3.4 (Piecewise affine and countably piecewise affine function). A function v : 0, — > 
M 2 is countably piecewise affine if v\ T is affine on every triangle T of a suitable locally finite 
triangulation of SI. If Q is a polygon and the triangulation is finite, then we say that v is 
(finitely) piecewise affine. 

In order to build the triangulation on which the function v of Theorem |1.1| is countably 
piecewise affine, we will use, on a subset of $7 of Lebesgue measure as close as we want to Jz? (f2), 
uniform triangulations into right triangles. The union of such triangles will be called a right 
polygon, according to the following definition. From now on, ei, e2 will be two fixed orthonormal 
vectors in M 2 . 
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Definition 3.5 (Right polygon and r-piecewise affine function). An open bounded set Cl' C R 2 
is called a right polygon of side-length r (or simply an r-polygon ) if it is a finite union of closed 
polygons whose sides are all parallel to e±, e2, and have lengths which are integer multiples of 
r > 0. Let now Cl' be an r-polygon, and consider a bi-Lipschitz function u : Cl' — > M?. For every 
side r C dQ.' , write T = U^LiTi where the T{ 's are essentially disjoint segments of length r. We 
say that u is r-piecewise affine on dCl' if for any such segment T and for any i, the function u 
is affine on Tj. 

Points in Cl will be denoted by z £ M 2 or by (x, y) £ R x R, with z = xe\ + ye2- We denote 
with B(z, r) the ball of center z and radius r and with T>(z, r) the square of center z, side length r 
and sides parallel to ei, e2- Moreover, the generic square of a collection of squares {V (z a , r a )} aG fq 
will be also sometimes denoted simply by T> a . Instead of working directly with triangulations, 
it will be convenient, in order to apply our method, to subdivide Cl into a countable and locally 
finite family of squares called tiling. 

Definition 3.6 (Tiling). Given an open, bounded set Cl, a tiling of Cl is a locally finite (in Cl) 
collection of closed squares {T> a (z a ,r Q )} a ^ whose union is contained between Cl and closCl and 
such that, Va ^ f3 € N, T> a V\T>^ is either empty, or a common vertex ofT> a and T>r, or a side of 
one of the two. Two squares of a tiling are said to be adjacent if their intersection is nonempty. 

Notice that a tiling of Cl can be either finite or countable and in particular it is surely 
countable if Cl is not a right polygon. 

It will be often useful to regard a given tiling of Cl as the union of the finite tiling corre- 
sponding to a right polygon Cl' CC Cl and a countable tiling of Cl \ CI' , locally finite in Q. Since 
these kinds of "sub-tilings" will be frequently used in the paper, we define them separately. 

Definition 3.7 (r-Tiling of a right polygon and tiling of (Cl,Q')). Given an r-polygon £1' , the 
r-tiling of $7' is the (unique) finite collection of closed squares {T>(z a , r )} a ^^f r ) whose union 
is equal to closil' and, Va 7^ f3 £ ^ (r), T> a n Pg is either empty, or a common vertex, or a 
common side ofD a and T>r- Given a bounded, open set Q and an r-polygon Q' CC Vl, a tiling 
of (fl,fl') is a tiling of CI whose restriction to CI' is the r-tiling of CI' . 

The 1-dimensional skeleton of a tiling will be called grid, according to the following defini- 
tion. 

Definition 3.8 (Grid). Let {T> a } a ^ be a tiling of CI. We call grid of the tiling the 1-dimensional 
set given by the union of the boundaries of the squares of the tiling. Each side (resp. vertex) of 
the squares of a tiling will be called side (resp. vertex,) of the grid. 

A possible definition of a piecewise affine approximation of u on a given r-polygon, which 
will be used in Section |4j is the following. 

Definition 3.9 ((CI', r)-interpolation of u). Let CI' CC Cl and {T> a } ae j? / r \ be an r -right polygon 
and its r-tiling. We call (Cl 1 , r)-interpolation of u the piecewise affine function v : Cl' — > v(Cl') C 
M 2 which coincides with u on the vertices of the r-tiling and, for each a G (r), is affine on the 



APPROXIMATION OF ORIENTATION-PRESERVING MAPS 



9 



two right triangles forming T> a and having as common hypothenuse the north-east/ south-west 
diagonal ofT> a . 

We conclude this section with a table collecting the main notation used in this paper. 



new 2 




a given open bounded set , 


V(z,r) 


square with center z, side length r 


u : f2 — 


■» A 


a given L bi-Lipschitz function , 




and sides parallel to ei , e 2 , 


i^(2x 


2) 


two by two real matrices , 


if 


Lebesgue measure on M 2 , 


|M| 




sup { \Mv\ : «| = l} , 




f-dimcnsional Hausdorff measure , 


i^(2x 


2;L) 


{M £^(2x2): DetAf > 0, 


hit A 


interior of a set ACM 2 , 






\M\ < L, M^ 1 < L] , 


clos A 


closure of A , 


ei, e 2 




two fixed positively oriented 


dA 


boundary of A C R 2 , 






orthonormal vectors in M 2 , 


n' cc n 


clos Q' C fi , 


B(z,r) 




ball with center z and radius r , 


d(A, B) 


inf{ \z - w : zeA,u)£iJ}. 



4. Approximation on the "Lebesgue squares" 
The aim of this section is to prove the following 

Proposition 4.1. For every e > there exists a right polygon Q e CC 0, of side length r such 
that the (Q £ ,r) -interpolation v : Q £ — > v(Q £ ) CR 2 is L + e bi-Lipschitz and satisfies 

A £ :=v(n £ ) CC A, (4.1) 



v 



u\\L°°(n £ ) + \\v~ - u~ ||l-(a £ ) + \\Du - Dv\\ LP (n e ) + \\Du~ - Dv ~ \\lp{A b ) < £ , (4.2) 

jS?(fi\n e )<e, ^(A\A e )<e, d{fl £ , R 2 \ ft) > 2r , (4.3) 
\/2r 

\\v-u\\ L oo (ne) < -^p- . (4.4) 

The reason why the piecewise affine interpolation of u will be injective on £l £ is that, for each 
square T> a of the r-tiling of £l £ , the function u will be uniformly close to an affine L bi-Lipschitz 
function on the nine squares around T> a . The linear part of each of these affine functions will 
be the differential of u at some Lebesgue point for Du inside T> a . For this reason, the squares 
of such r -tiling will be called "Lebesgue squares". 



The plan of this section is the following. Section 4.1 contains Lemma 4.2, which is the main 
ingredient in the proof of Proposition |4.1| Indeed, Lemma 4.2 says that, when on a square Du 
is close in average to an L bi-Lipschitz matrix M, then u is close in L°° to an affine function um 
with Dum = M. Then, in Section |4.2[ we will determine Cl £ as a suitable union of squares of 
an r-tiling on which Lemma |4.2| holds and provides a sufficiently strong L°° estimate. Finally, 



in Section 4.3 we show that the (fi e , r)-interpolation of u satisfies the required properties. 



4.1. An L°° Lemma. We are now ready to begin the proof of Proposition |4,1[ starting from 
the following fundamental lemma. Here and in the following, by ^# (2 x 2; L) we denote the 
set of the two by two invertible matrices such that the affine map z 1— > Mz is L bi-Lipschitz. 
Moreover, f2 and u will always be a set and a function as in the assumptions of Theorem |1.1| 
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Lemma 4.2. For any rj > there exists 5 = 5(n) > such that, if z G Q, M G ^#(2 x 2; L) 
and p > are so that D(z,p) CC 0, and 

I \Du{z)-M\dz<5, (4.5) 

JV(z,p) 

then there exists an affine function um '■ K 2 — > M 2 with Dum = M and such that 

\u{z)-u M {z)\<np Vz£V(z,p). (4.6) 

Proof. Up to a translation, we are allowed to assume for simplicity that z = u(z) = (0,0) G M 2 . 
Let us then call, for a big constant R to be specified later, 

rp/2 



B 1 := jx G [-p/2, p/2] : J* \Du(x, t) - M\ dt < pR5^ , 
B 2 := <|y G [ — p/2, p/2] : ^ |.Du(t, y) - M\ dt < pRS^j . 



-p/2 

Notice that, since u is bi-Lipschitz on $7, then so are its restrictions to the horizontal and vertical 



segments of the square T>(0, p). Hence, the above integrals make sense for every x and y. By (4.5) 
and Fubini-Tonelli Theorem, we readily obtain 

JZ> 1 ([- p/2, p/2] \ B 1 ) < | , Jf 1 ( [ - p/2, p/2] \ B 2 ) < | . (4.7) 

Define now um{z) = Mz, and <p(z) = u(z) — um{z). For any x\, X2 G B 1 and y%, yi G -B 2 we 
immediately get 

|^(^i,yi) - <p(»2,y2)| < |<p(xi,?/i) - cp(x 2 ,yi)| + |cp(x 2 ,yi) - ^(^2, 2/2) | 

(4.8) 



< \Du(t,yt) - M\dt+ / |D«(x 2 ,t) - M| dt < 2pR5 . 

J xi J yi 



Let now (x, y) G T>(z, p) be a generic point. By (4.7), there exist x\ G B 1 and y\ G -B 2 so that 

\x-xi\<—, \y - yi\ < — , 

and since u and um are L bi-Lipschitz, thus <p is 2L-Lipschitz, we get 

\(p(x,y)-(p(x liyi )\ < 2 lMi. (4.9) 



Let finally (x,y) and be two generic points in D(z,r). Putting together (4.8) and (4.9) we 

immediately get 

\v{x,y) - f(x,y)\ < + 2 P R5 ^ VP, 

where the last inequality is true up to take R big enough and then 5 small enough. Since 
(p(0, 0) = 0, this concludes the proof. □ 
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4.2. A large right polygon made of Lebesgue squares. In this section we show that, 
for any 77 > 0, it is possible to construct a right polygon Q„ CC of side length r v such that 
Jz? (Q\£l r) ) < 77 and such that, for any square V(z,r v ) of the r^-tiling of O^, the assumption (4.5) 
of Lemma 4.2 holds on the bigger square T>(z, 3r v ). As we will show in Section 4.3 , if we choose rj 
and then r v small enough, the corresponding r J) )-interpolation of u satisfies the requirements 
of Proposition 4.1 Then, Q.^ will turn out to be the right polygon of Lebesgue squares we are 
looking for. The goal of this section is to show the following estimate. 

Lemma 4.3. For every 77 > there exists a constant r = r{rf) > and an r-polygon fi^ CC 
such that J§? < rj and each square of the r -tiling {T>(z a ,r)} r r \ satisfies the following 

properties, 



V(z a ,3r) CCn \/a£J?{r) 

f \Du(z) -M\dz< 5(rj) for some M = M(a) G Ji (2 x 2; L) . 



(4.10) 
(4.11) 



>V(z a ,3r) 

Proof. We start selecting some ro = 7*0(77) > and an ro-polygon f!o CC H such that Jz? (O \ 
^0) < 77/2 and each square of the ro-tiling of f^o satisfies (4.10). Then, for every r such that 
ro G rN, we can regard Qq also as an r-polygon, and consequently call {V(z a , r )} a6 / ( r ) its 
r-tiling. We define the set 



J (r) :-- 



a^J (r) : f \Du - M\ < 5 for some M = M{a) ei'(2x 2;L)|, 



where 5 = 5{rj) is given by Lemma 4.2 and we let 



n 



v •" 



Since property (4.11) is true by construction, to conclude the proof it is enough to select a 
suitable r = r{rj) in such a way that if (O \ 0„) < 77/2. 

To do so, we apply the Lebesgue Differentiation Theorem to the map Du finding that, for 
J2?-a.e. z G $7o> there exists r(z) > such that V(z,4r(z)) C Q and 

5 



V{z,p) 



\Du(w) — Du(z) \ dw < 



V0 < p < 4r(z) 



We can then choose r = r(rj) so small that the set A(r) : = {z £ Slo : r{z) < r} satisfies 

&(A{r)) < 77/2. (4.12) 
We now claim that, for each a G J^o( r )i 

P(z a ,r)^i(r) =0 a£j(r). (4.13) 
Indeed, letting M = Du(z) for some 2 G P(^ a ,r) \ A(r), by definition of A{r) and r(z) we get 



T>(z a ,3r) 



\Du-M\ 



1 

972 



\Du-M\ < 



T>(z a ,3r) 



9r 2 



16 f 

\Du-M\ = —f \Du-M\ 
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thus (4.13) is obtained. As a consequence, by (4.12) we have that 

££ ( Q \ n v ) = ££ ( I J 
and, as we noticed above, this concludes the proof. 



V(z a ,r) \ <J?(A{r)) < 



□ 



4.3. Affine approximation of u on Lebesgue squares. In this section we complete the 
proof of Proposition |4.1| At this point the proof reduces to show that, provided we choose 77 



small enough, the (f^, r)-interpolation of u on the right polygon as in Lemma 4.3 satisfies 



the properties of Proposition 4.1 



Proof of Proposition J^.l; Let e > be a given constant. Then, let 77 = 77(e) be a sufficiently 
small constant, whose value will be precised later. Define now 5 = (5(r/(e)) as in Lemma 4.2, and 
define also r = r (77(e)) and f2 e = fi^M according to Lemma 4.3 We will show that the right 



polygon ft £ fulfills all the requirements of the proposition as soon as 77(e) is small enough. To this 
aim we call, as in the statement, v : il £ 
on the right polygon Q E . 



► A £ the (0 £ , r)-interpolation of u (see Definition 3.9) 
Let us briefly fix some notation which will be used through the proof. For any a 6 ^(r), 



we define M a £ ^ (2 x 2; L) so that (4.11 ) holds. Applying Lemma 4.2 with p = 3r, we get an 
affine function u a : M 2 — > M 2 such that Du a = M a and 



\u — u a \ < 3777" on V(z a , 3r) 
Figure [4] depicts the functions u, v and u a . 



(4.14) 



D 



C 



A 



B 




u a {A) 



u a {B) 



Figure 4. The functions u, v and u a on a square. 



We can then start the proof, which will be divided in some steps for clarity. 

Step I. For any a G (r), v(T>(z a , r)) C u(V(z a ,3r)) . 

Take a € (r). Keeping in mind ( 4.14| ) and recalling the definition of v, we get that 



v(V(z a ,r)) C B(u a (T>(z a ,r)),3rjrJ , 



(4.15) 
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where for any set I C I 2 we denote by B(X, r) the r-neighborhood of X. Similarly, we get that 

x : B(x, 3nr) C u a (V(z a , 3r)) | . 

Hence, the step is concluded if 

B(u a (V(z a ,r)),6rir\ C u a (T>(z a ,3r)) , 

which in turn, recalling that Du a = M a £ j% (2 x 2; L), is true as soon as r\ < (6L) _1 . 

Observe that, as an immediate consequence of this step and (4.10), we have A £ CC A, that 



is, (4.1) holds. 

Step II. Injectivity of v. 

Take a 6 <# {r). Applying again (4.14) as in Step I, we deduce that v is injective on 
T>(z a , 3r)nf2 e as soon as rj < (6L) _1 . To conclude that v is injective, then, we have to show that 
v(V(z a , r))nv(T>(zf3, r)j = if T>(z a ,r) and T>(zp,r) are two non-adjacent squares of the tiling of 
O e . And in fact, if T> a and T>p are non-adjacent, then the fact that v(T>(z a , r)) r\v(T>(zp, r)) = 
for a small t] follows as an immediate consequence of (4.14) and ( |4. 15 ) arguing as in Step I. 

Step III. Estimate on \\v — u\\uoo(qa and on Wv^ 1 



u IU°°(A e )- 



Fix a generic square T> a of the r-tiling of O e , and observe that \\u a — u\\Loo(p\ < 3r]r 
by (4.14). Moreover, v and u a are both affine on each of the two right triangles on which T> a is 
divided, and since on the vertices of these triangles v equals u, again by (4.14) we deduce also 
II v ~ u a\\L oa {v a ) ^ 3nr. Thanks to these two estimates, we deduce 



\v ~ w|U<»(fi e ) = sup \\v 

ae/(r) 



U>\\L°°(V a ) < SU P \\ v 



U-aWl^iVa) + ll^a — u \\L°°(V a ) 



(4.16) 



< 6r/r < 



4L' 



where the last inequality is true as soon as rj, hence also r, is small enough. 

Since we have already proven that v is injective, the L°° estimate for the inverse maps is 
now a simple consequence. Indeed, taking a generic point w = v(z) £ A e , with z E Q £ , by (4.16) 
we have 



so that 



(w) — v 1 (w)\ = \u 1 (v(z)) — u 1 (u(z))\ < L\v(z) — u(z)\ < 



-ill <- 
\\l°°(a £ ) S 



(4.17) 



Step IV. Estimate on \\Dv — Du\\lp(q\. 

Let us start observing that, since by construction \Du\ < L and \Dv\ < \/2L, one has 



\Dv - Du\\ p LP(Ue) 



E ii 

p-1 



< 3L 



< (3L 



Dv - Du\\ p LP{Va) 



E W Dv ~ Du \\L\V a ) 



(4.18) 



p-i 



E 



\Dv - Du a \\ L i(v a ) + \\ D 



Du\ 



LHV a ) 
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By (4.11 ), we already know that for each a G ,f (r) 



Du- Du ol \\ L i (Va) = f \Du- Du a \ <9r 2 f \Du - M Q \ < 95r 2 = 9S\V a \ . (4.19) 

JV(z a ,r) Jv(z a ,3r) 



(z a ,3r) 



Let us then concentrate on \\Dv — Du a \\£i(p a y Consider the triangle T = z±z 2 z 3 , being 

z 1 = z a +(- r/2, -r/2) , z 2 = z a + (r/2, -r/2) , z 3 = z a + (r/2, r/2) . 

Since both t> and u a are affine on T, then in particular — Du a is a constant linear function 
on T. Recalling again ( 4.14| ), let us then calculate 

{Dv\ T - Du a ) (rci) = («(z 2 ) - - {u a (z 2 ) - u a {zi)) 

= (u(z 2 ) - u(zi)) - (u a (z 2 ) - u a (zx)) 

and similarly 

(Dv\ T - Du a )(re 2 ) = (v(z 3 ) - v (z 2 )) - (u a (z 3 ) - u Q (z 2 ) 

We deduce that \\Dv — Du a \\ ^00^ < Q^/2r]. We can argue in the same way for all the different 
triangles in which V(z a , 3r) D Q £ is divided, thus getting 

\\Dv- Du a \\ Lao{v{zai3r)nne) < 6V2t] < 977. 

Inserting this estimate and ( |4.19[ ) into ( |4.18[ ), we get 



< 6r]r , 



< 6r/r . 



(4.20) 



\Dv - Du\\ p LP(Qs) 



^{SLf-'d^ + v) J] \V a \ = (3Ly-h[ V + 6)\n E \<(~Y , (4.21) 



where again the last inequality holds true as soon as rj, hence also 5, is small enough. 
Step V. Bi-Lipschitz estimate for v. 



Take a point z G T>(z a , 3r) n Recalling that u a is L bi-Lipschitz and (4.20), we get 



1 
L 



9t] < \Dv(z)\ <L + 9r/. 



(4.22) 



Let then z, z' G H £ be two generic points, and assume that z G V(z a , r). If one has z' G V(z a , 3r), 
then an immediate geometric argument using the definition of v and (4.22) yields 



1 



- -9?? < |w(z) < {L + 9r])\ z — z . 



(4.23) 



On the other hand, assume that z' ^ T>(z a ,3r), so that \z — z'\ > r. In this case, the L°° 
estimate (4.16) gives 

\v{z) - v(z')\ < \u(z)-u(z')\ + \v(z)-u(z)\ + \v(z) -u(z')\ < (L + 12r))\z - z'\ , (4.24) 
and similarly 

\v(z) - v{z')\ > \u(z) - u(z')\ - \v(z) - u(z)\ - \v(z') - u(z')\ > (- - 1277 J \z - z'\ . (4.25) 



Putting together (4.23), (4.24) and (4.25), provided that r\ is small enough we conclude that v 
is L + e bi-Lipschitz. 

Step VI. Estimate on \\Dv~ 1 — Du~ 1 \\lp(A e )- 
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Fix a generic a £ ^ (r), and recall the elementary fact that, given two invertible matrices 
A and B, one always has |£> _1 — A~ 1 \ < | \B — A\. Since by construction u and u a 
are L bi-Lipschitz, and Du a is constant on T> a , then Step I, (4.11) and (4.19) ensure that 



\Du 



Du <x \W{v(D a )) 



\Du- l (z) - Du^ 1 (z)\dz 



u(D(z a ,3r)) 



v(D(z a ,r)) 

D«(tt -1 (z)) -M a \dz < L 4 



\Du(w) — Mai dw 



T>(z a ,3r) 



9r 2 L 4 



/ 



I Du — M a I < 9r 2 L 4 5 = 9L 4 5 \V a 



>V(z a ,3r) 

On the other hand, again using \B~ 1 — A~ 1 \ < 1 \B~ l \\B — A\, the fact that Du a is constant 
on D a , the fact that u a is L bi-Lipschitz by definition while v is (L + e) bi-Lipschitz by Step V, 
and (4.20), we readily obtain 

,-i 



\Dv~ 



Du 



-ill 

a \\L°°(v(V a )) 



< L(L + e)9ri < 18L 2 r/ . 



We can then repeat the same argument as in (4.18) to get 

\P~ 1 



\Dv~ 



Du-% p[Ae) <(3L) 



E ii^" 1 



<(3L) 



p-i 



Du a \\v(v(V a )) 



\Du- 



Du 



\L^(v(V a )) 

(4.26) 



18^7/1 A e | +9L 4 5|O e |J < 

where as usual the last estimate holds possibly decreasing r] and then also 5. 
Step VII. Conclusion. 



Let us now conclude the proof of Proposition 4.1 by checking that f2 e fulfills all the require- 
ments of the statement. The fact that v is L + e bi-Lipschitz is given by Step V. The validity 
of fl3~I] ) has been observed in Step I. The estimate (J42j) just follows adding fl4.16[ ), \A.17\ , ( |4.21[ ) 
and (4.26). Concerning (4.3), the facts that f2 \ J7 e is small and that cf(f2 £ ,IR 2 \ O) > 2r are 
given by Lemma 4.3, while the fact that also A \ A e is small is immediate by the bi-Lipschitz 
property of u and the L°° estimate (4.6) of Lemma 4.2 Finally, (4.4) is immediate provided 
that we choose rj < \/2/(36L 3 ), since (4.16) ensures that \\v — fct||L°°(fi s ) ^ 6r/r. □ 



5. Approximation out of "Lebesgue squares" 

In this section we complete the proof of Theorem |1.1[ defining the countably piecewise 
affine approximation of u out of the large right polygon f2 e of "Lebesgue squares" constructed 
in Proposition 4.1 Following the scheme outlined in Section [2j the construction is carried out 
in three steps: the covering of \ Q £ with a suitable (locally finite) tiling, the definition of a bi- 
Lipschitz piecewise affine approximation of u on the grid of the tiling and, finally, the extension 



of the approximating function to the interior of the grid by means of Theorem 2.2 The main 
result of this section is the following. 
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Proposition 5.1. Let v £ : f2 £ 

tion 



A E be a piecewise affine bi-Lipschitz function as in Proposi- 



ti 



Then, there exists aC\L^ bi-Lipschitz countably piecewise affine function v £ : fi\fl s — > 
A \ A £) where C\ is a geometric constant, such that v £ = u on dVl and v £ = v £ on dfl £ . 

We can immediately notice that Theorem |1.1| follows as an easy consequence of Proposi- 
tions Hj] and ED 



Proof of Theorem \l.l\ Take e > 0, and apply Proposition 4.1 to get an r-polygon S7 £ CC £1 
and a piecewise affine bi-Lipschitz function v £ : £1 £ — > A e . By Proposition 5.1, we have a C\L 4 
bi-Lipschitz function v £ : O \ f2 e — > A \ A e , so we can define the function v : — > A as v = v £ 
on £l £ and v = v E on O \ f2 e . Since v £ (resp., v £ ) is bi-Lipschitz with constant L + e (resp., 
C\L ), and v £ = v £ on <9£\, we have that v is a bi-Lipschitz homeomorphism with constant 
C\L . Moreover, by construction v is countably piecewise affine and it is orientation-preserving, 
since so is u and v = u on d£l. We are then left with showing that v satisfies (1.1), and by (4.2) 
this basically reduces to consider what happens in f2 \ fi e . Since v £ is bi-Lipschitz with constant 
CiL 4 , by (4.3) we clearly have 



\\Dv 
and similarly 



Du 



\ L p(n\n e ) < \\Dv-Du\\ Loo{n \ nE) \n\n £ \ 1/p ^(L + dL^e 1 ^ 
\\Dv- 1 - Du' 1 || LP(A \A £ ) <{L + C 1 L 4 )e 1 / p . 



(5.1) 



(5.2) 



Concerning the L°° estimates, since O \ flJ < e then for every z G $7 \ Q, £ there exist z' 6 Q £ 
such that \z — z'\ < y/e/n, thus by (4.2) we find 

\v(z) — u(z)\ < \v(z) — v(z')\ + \v(z ) — u(z')l + l u 0O ~~ 

< (L + dL 4 ) ^1 + [|« 6 - u|U» ( n e ) < (i + Ci^ 4 ) a/J + e • 
Arguing in the same way to bound \v (w) — -u _1 (u>)| for a generic toe A\A £ yields 

lb - «lll/~(n\n e ) <(■£- + Ci^ 4 ) W^; + e , H^ 1 - ||i,oo(A\A e: ) < + Ci-^ 4 ) + ^ • (5-3) 

Putting together ( |5.1[ ), (5.2) and (5.3), we find that t> satisfies ( |1.1[ ) as soon as e is chosen 
small enough, depending on e. Hence, we have found the countably piecewise affine approxima- 
tion as required. Concerning the smooth approximation, its existence directly follows applying 
Theorem 



2.1 



thus we have in particular C2 = 706^ 



7/3 



□ 



We have now to prove Proposition 5.1 To do so, let us fix some notation. Recall that Q £ is an 
r-polygon for some r = r(e). We will then start by selecting a suitable tiling {Vj = V(zj,rj)}j £ fq 
of (O, O e ), according with Definition 3.7 This means that {T>j} is a tiling of SI whose restriction 
to £l £ coincides with the r-tiling of f2 £ . The only requirements that we ask to {T>j} are the 
following, 



Vj : ciosP,- ndn £ / 



PjCcO VjeN. 



(5.4) 
(5.5) 
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Notice that (5.4) is possible thanks to (4.3), while (5.5) basically means that the tiling has to be 



countable instead of finite, and the squares have to become smaller and smaller when approaching 
the boundary of Q. Of course, in the particular case when f2 itself is a right polygon, instead 



of (5.5) one could have asked the tiling to be finite (we will discuss this possibility more in detail 



in Remark 5.7). 



Since it is of course possible to find a tiling of (f2,O e ) which satisfies (5.4) and (5.5), from 



now on we fix such a tiling, and we denote by Q its associated 1-dimensional grid according to 
Definition 3.8 Moreover, we set Q' = Q n (f2 \ closf2 e ), which is the part of the grid on which 
we really need to work. Notice that Q' is a 1-dimensional set, made by all the sides of the grid 
Q which lie in O \ clos il e . 

Let us call now w a the generic vertex of Q! , hence, the generic vertex of the grid Q which 
does not belong to Q £ (however, notice that w a may belong to dQ £ , since by definition Q' is 
relatively open around dQ £ ). Each vertex w a is of the form w a = Zj + (±rj/2, ±rj/2) for some 
j, and it is one extreme of either three, or four sides of Q. To shorten the notation, we will 
denote the other extremes of these sides by w l a with 1 < i < i(a), being then i(a) £ {3, 4}. 
Finally, we will denote by t a the minimum of the lengths of the sides w a w l a . Observe that if 
w a ^ dQ £ , then w a is one extreme of either three or four sides of Q' C Q. On the other hand, 



if w a G <9S\, then by (5.4) it is one extreme of four sides of Q, either one or two of these four 
sides lies in Q', and i a = r. 

" we 



Thanks to Theorem 2.2 to obtain the piecewise affine function v £ of Proposition 5.1 



essntially have to define it, in a suitable way, on the 1-dimensional grid Q' . To do so, our main 



ingredients are the following two lemmas. The first one (Lemma 5.5) states that, on any given 
segment of O, u can be approximated as well as desired in L°° with suitable piecewise affine 
3L bi-Lipschitz functions. This is of course of primary importance to define the piecewise affine 
approximation v £ of u on the sides of the grid Q', but it is still not enough. In fact, we have to 
take some additional care to treat the "crosses" of Q' (that is, the regions around the vertices), 
in order to be sure that our affine v £ on Q' remains injective. This will be obtained thanks to 
the second Lemma 15.61 

To state the next two lemmas, it will be useful to introduce some piece of notation. 



Definition 5.2 (Interpolation of u). Given a segment pq CC 0, let {-Zi^j+i}o<j<A r ^ e N essen- 
tially disjoint segments whose union is pq, with zq = p and zn = q. For any such subdivision 
of the segment, will call interpolation of u the finitely piecewise affine function u pq : pq — > M? 
such that, for any < i < N — 1 and any < t < 1, 



u 



Zi + t{z i+ i - z^ = u(zi) + t(u(z i+ i) - u(Zi)) . 



Definition 5.3 (Adjusted function and crosses). Let {^a}a&i be a sequence such that for any 
a one has 3L£ a < £ a . For any a G N and any 1 < i < i(a), we define as the biggest number 
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such that 



u(w a ) - u(w a + f a « - w a )) 



U[W n 



< £a if WaWa C Q' , 



We w«ZZ call adjusted function the function u a dj : Q — >■ M 2 defined as follows. First of all, we 
set u a dj = v £ on Q\Q' . Then, let w a wp be a side of Q' , thus being wp = w l a and w a = wp for 
two suitable i and j . We define 



u&djiwa + tiwp - w a )) := 



U(w a ) + g^(w a + C a (wp - W a )) - u(w a ) 
u(w a + t(wp - w a )) 

u{wp) + (u(w/3 + ^{Wa - Wp)) - U(wp 



in (0,£J, 
in (1-&1). 



/n words, for any side in Q! , u a dj coincides with u in the internal part of the side, while the two 
parts closest to the vertices w a and wp are replaced with segments. Moreover, for any vertex w a 
of Q! we will define its associated cross as 



i(a) 

1=1 



: < t < 



Remark 5.4. Some remarks are now in order. First of all, since u is L bi-Lipschitz on f2, and 
also v e is L bi-Lipschitz on any segment w a w a C Q \ Q' , by the choice 3L£ a < l a we directly 
deduce that < < 1/3 for any a and any i. Thus, two different crosses have always empty 
intersection. For the same reason, each of the i{a) extremes of the cross Z a has a distance at 
least £, a /L from w a . Finally, for all different a and (3 one has B(u(w a ),^ a ) r\B(u(wp),^p) = 0. 
Indeed, assuming without loss of generality that £ a > ip, we have \u(wp) — u(w a )\ > £a/L. And 
as a consequence, £ a + £p < t a j (3L) + £p/(3L) < 2£ a / (3L) < \u{wp) — u{w a )\. 



Lemma 5.5. For every segment pq CC f2 and every 5 > 0, there exists a function u pq : pq 
which is a 4L bi-Lipschitz interpolation of u with the property that \\Up q — ■u||l°o(p 9 ) < S. 



A 



Lemma 5.6. There exists a sequence {£, a }aeN such that the associated adjusted function u a dj '■ 
Q — > R 2 is 18L bi-Lipschitz and w a dj(Q) Q A. 



Before giving the proof of Lemmas 5.5 and 5.6, we show how they enter into the proof of 



Proposition 5.1 



Proof of Proposition 5.1; To define the searched function v £ : Q, \ Q 
» M 2 be an adjusted function according with Lemma 



5.6 



— > A \ A e , let first 
corresponding to the 



^adj : Q — 

sequence Our strategy will be first to define a suitable piecewise affine and injective 

function u'^- : Q — > A, coinciding with n a( jj near the vertices w a , and then to obtain v £ 
extending u'^ in the interior of each square making use of Theorem 2.2 We divide the proof 
in some steps. 



Step I. Definition of u' ad - : Q — > A. 
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First of all, we define it ad j = u a ^ = v £ on Q\ Q' . Then, we consider a generic side w a wp C Q' 
and define pq the internal segment of the side w a wp, that is, p and q are the extremes of the 
segment w a wp \ (Z a U Z$). Taking now a small constant 5 = 5(a, /?), to be fixed later, we set 
u adj = ""adj on w a W/3 n (Z a U Zp) , and u' ad - = u pq on pq, where u pq is given by Lemma 



5.5 



By definition, it is clear that u' ad j is a continuous, countably piecewise affine function on 
Q. Moreover, since the different constants 5(a, f3) can be chosen arbitrarily small, each one 
independently of the others, and any internal segment pq is compactly supported in f2, one can 
clearly assume that u' ad -(Q) C A. In addition, since u' ad j is obtained gluing the 4L bi-Lipschitz 
functions u pq and the 18L bi-Lipschitz function n a( jj, we clearly have that u' ad - is 18\/2-^-Lipschitz 
(but, a priori, not bi-Lipschitz and not even injective!). To conclude the proof, we will then 
show that in fact u ad j is bi-Lipschitz (hence, in particular, injective), and eventually we will 
extend ii ad j to the interior of the squares of the tiling (hence, to the whole Q \ Q £ ). 

Let us then fix two points z, z' £ Q. In the next Steps ILTV we will show that, provided 
that the constants S(a, f3) are chosen small enough, 

|<dj(*)-<ij(2')|> ^f\ z -A- (5-6) 
We consider separately the different possible reciprocal positions of z and z' . 

Step II. The case in which z £ pq C w a W/3, z' £ w a wp . 

In this observed in Remark |5.4[ we know that \z — z'\ > £, a /L. Thus, there are two 

subcases. If z 1 does not belong to any internal segment (hence, either z' belongs to some cross, 



or z' G O e ), then u' &A ^{z') = Ug^j(z') and then by Lemma 5.6, provided that we choose 



r / ns . min{£ Q , £ } 

jjg£2 > ( 5 - 7 ) 



we have 



Kdj( z ) - U U]{ Z ')\ = |«Jff(*0 ~ «adj(^)| > hadj(^) - «adj (z') \ ~ \u & pq {z) - U ad j(z)\ 

= |«adj(z) - u adj (z')\ - \u 5 pq (z) - u(z)\ > \z - z'\ - 5 (a, /3) 

> \Z — Z \ ^7T > 



18L 1 1 36L 2 ~ 36L 



so that (5.6) is proved. 

Consider now the other subcase, namely, when z' belongs to some other internal segment 
p'q' C w a iwpi. In that case, since by construction and (5.7) it is \z — z'\ > 36L5(a, [3) and 
\z — z'\ > 36L5(a' , /?'), one directly has 



> 1 |* - z'\ - 5(a, f3) - 5{a', tf) > ^\z - z'\ 



hence again (5.6) is established. 

Step III. The case in which z G pq C w a wp, z' 6 w a wp . 

The second case is when z still belongs to an internal segment contained in the side 
w a wp C Q', and also 2' belongs to the side w a wp. In particular, if also z' is in the internal 
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segment pq then we already know the validity of (5.6) because w' ad j(z) = u pq {z) and u' ad -(z') = 
UpJz'), while Up q is 4L bi-Lipschitz. Therefore, we can directly assume that z' S w a p, being the 
case z' G qwp clearly the same. 

By Definition 5.3 we know that it ad j(z') = Wadj^') lies in the segment u(w a )u(p), which is a 
radius of the ball B(u(w a ), £ Q ) . Hence, for any point s outside the same ball, a trivial geometric 
argument tells us that 

- u(p) \ + \u{p) - u adj (z')\ 



^adjl 



> 



(5.8) 



Notice now that it is not true, in general, that u ad j (z 



u, 



pq 



However, recalling that u pq is an interpolation of u, by Definition 
segment whose both extremes are out of the ball. Thus, if 
decreasing <5(a, /3) it is surely true that 



z) lies outside the ball B(u{w a ), £ a ) . 

5 ( 

pq\ 



5.2 



u 



(w a )\ < £ a , up to possibly 



£a - \u pq (z) - u(w a )\ < \u pq (z) - U(p) \ . 

Putting this observation together with (|5.8|) we readily obtain that 



K 9 0) -^adjOOl > 



\u s p Jz) - u(p) \ + \u(p) - M adj (2: / )| 



U pq( Z ) ~ Up q {p) \ + \u ad j(p) ~ U ad j(z')\ 



recalling that n ad j(p) = u s pq (p) = u{p) (of course, by selecting 5(a, /3) small enough, we could 
have used any number greater than 3, instead of 4, in the above estimate). Therefore, since Up q 
is 4L bi-Lipschitz while n ad j is 18L bi-Lipschitz, we readily obtain 



*adj 



4 adj 



CO I 



> 



\ U pq( Z ) 

\z — p\ 
16L 



+ 



«adj(^')| > 
\p - 



72L 



> 



\ U pq( Z ) 



72L 



uL{p)\ 



|«adj(p) - ^adj(2')| 

~~ 4 



recalling that z, p and z' are aligned. Hence, (5.6) is checked once again also in this case. 
Step IV. The case in which neither z nor z' are in some internal segment. 



Thanks to Step II and Step III and by the symmetry of the inequality (5.6), we are left 
to consider only the situation where no one between z and z' is inside some internal segment. 
In other words, both z and z' must be either in Q\ Q' or in some cross. By the definition 
of it ad j, this means that « ad j(z) = "u ad j(z) and u' ad ^(z') = n ad j(z'). And thus, since -u ad j is 18L 
bi-Lipschitz thanks to Lemma 5.6, the validity of (5.6) is already known. Summarizing, we have 
shown the validity of (5.6) in any possible case, and this means that the function u'^ : Q — > A 
is injective and 72L bi-Lipschitz. 

Step V. Conclusion. 

We have now to define the piecewise affine and bi-Lipschitz function v £ : £1 \ 0, £ — > A \ A e , 
matching u on <9Q and matching v e on dVt e . To do so, consider each square T>j of the tiling 
contained in O \ f2 e . The function n ad j is 72L bi-Lipschitz from dT>j to a subset of A, then 
by Theorem 2.2 it can be continuously extended to a piecewise affine bi-Lipschitz function of 
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the whole square T>j, with bi-Lipschitz constant C^J2 4: L . Define v £ as the countably piecewise 
affine function on Q \ £l e which gathers all these extensions on all the squares Dj C Q \ £ of 
the tiling. 

For each square Vj C Q \ il e , we clearly have that d(y e (T>j)) = u'^idVj). This yields that 
v e is injective. Moreover, by continuity it is clear that v e = u on d£l, and by construction v e = v e 
on d£l £ . As a consequence, v £ : f2 \ £l £ — > A \ A e fulfills all our requirements. In particular, 
one has C x = 72 4 C 3 . □ 

Let us now make a simple observation, which will be useful in the sequel. 
Remark 5.7. Assume that $7 is a right polygon of side-length f and that u is f -piecewise affine 



on dQ, according to Definition\3.5[ Then, consider the r -polygon f] £ given by Proposition 4-1 



By the construction of Section it is not restrictive to assume that r £ rN, and that f2 e is a 



subset of the r -tiling o/O. Therefore, we can repeat verbatim the construction of Proposition 5.1 



using, as tiling, the r -tiling of Q. Notice that in this case assumption (5.5) is not valid -see 



the remark right after (5.5)- but in fact if VL is an r-polygon, and u is r-piecewise affine on 



d£l, there is no need for the tiling to use smaller and smaller squares at the boundary. As a 



consequence, the bi-Lipschitz approximation v e provided by Proposition 5.1 is (finitely) piecewise 
affine instead of countably piecewise affine. Observe that the assumption that u is f -piecewise 
affine on dtt is essential, because otherwise the map v e would not coincide with u on dQ. 



To conclude the proof of Theorem we then only need to give the proofs of Lemma 5.5 
and of Lemma 15.61 



Proof of Lemma \5. 5\ Let p > be a small number, to be fixed later. Define then to = 0, zq = p 
and then recursively 

t i+ i := max jl >t>U: \u(zi) - u(p + t(q - p)) \ < pj , z i+ i :=p + t i+1 (q - p) . 
In this way, we have selected a finite sequence of points zq = p, z\, . . . zm = q in the segment 

'pq 



pq, where N = N(p,q,p). We can then already define the function u & pq by setting, for any 



< i < N - 1 and any < t < 1, 



u 5 
u pq 



\Zi + t(z i+ l - Zi)^j = u(zi) + t(u(zi + i) - u(zij) . 



5.2 



Hence, u 5 pq is the interpolation of u associated with the points {zi}, according with Definition 
The function u pq is by construction finitely piecewise affine and L-Lipschitz. By the uniform 
continuity of u in pq it is also clear that the bound \\u — ^H^oo^ < 5 holds true as soon as p 
is small enough. To conclude, we have thus only to check that 

K ? W-< 5 (z')|>^-^'| (5.9) 

for all z, z' in pq. If both z and z 1 belong to the same segment z%Zi+\, then the estimate is true 
because u pq is affine on that segment and u is L bi-Lipschitz. 

Assume then that z S ZiZi+i and z' 6 ZjZj + \ with j > i. lij = i+1, that is, z and z' belong to 
two consecutive segments, then by the definition of the points Zi the angle u pq {z)u 5 pq {zi + i)u pq (z') 
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is at least tt/3. Hence, 
\utlz)-ut Q {z!)\>- 



Upq(z) ~ U S pq (z i+1 )\ \u 5 pq (z i+1 ) - U 5 pq ( Z ') 



+ 



\u(zi) - u(z i+1 )\ 



2\zi 



+ \Z, 



\u(z i+ i) - u(z i+2 ) 



i+1 



2 z 



> 



i+1 — Zi + 2\ 



\z — z 
2L 



so that (5.9) is checked. 

Instead, let us see what happens if j > i + 1. In this case, since u pq (z) G u(zi)u(zi+i) and 
for alH > i + 1 one has u(zi) ^ B(u(zi),p) U B(u(zi + \), p), an immediate geometric argument 
ensures that \u pq (z) — u pq (z')\ > \f3pj2. As a consequence, we have 



\u(zi) - u(z j+ i) \ < \u pq (z) - 
which yields 



"pq< 



z')\+2p<(\ + ^ y/£j \u s pq (z)-u s pq (z')\ < 4\u s pq (z)-u 5 pq (z')\ , 

u{Zi) - u(Zj +1 )\ |, 



Zj+1\ 



4L 



\z — z 
4L 



hence (5.9) holds true also in this case and we conclude the proof. 



□ 



Proof of Lemma \5. 6} Let us take a vertex w a of the grid Q' . Take then a constant £ a < £ a /(3L), 
with £ Q = £ a /(3L) = r/(3L) if w a E <9£\, while if w a ^ the inequality can be strict. In 
particular, it is admissible to ask that for any a one has 

r 



2L 



(5.10) 



Define now ^ as in Definition 
w a 6fi \ <9f2 e , then we have 



5.3 



and, for any 1 < i < i(a), let Pi = w a + € a { w a ~ w a)- If 



u{w a )u{pi) CC A VI < i < i(a) , 



(5.11) 



up to possibly decrease the value of £ Q . Instead, if w a G 9J7 e , then (5.11) is already ensured 



by (4.4) and (4.3) in Proposition 4.1, without any need of changing £ a . 



We introduce then the adjusted function tt a( jj of Definition 5.3 to obtain the thesis, we need 
to check that it fulfills the requirements of Lemma 5.6 Thanks to (5.11 ), we already know that 
Undj ■ Q — > A. Hence, all we have to do is to check that 



18 L 



< |tiadj(z) - ^adjOOl < l8L\z - z'\ . 



(5.12) 



for all z, z' E Q. We will do it in some steps. 



Step I. For all a, yclos B[u(w a ), £ Q )J = Z a . 

We start observing an important property, that is, for any a and for any z £ Q we have 
that |n a dj(^) — ii(u> Q )| < £ a if and only if z S Z a . In fact, if z € Z a then z G w a pi for some 
1 < i < i(a), and since u a dj is affine in the segment w a pi, while |n a dj(pi) — u(w a )\ = Co, then of 
course |n adj (z) - u(w a )\ < 

On the other hand, assume that z ^ Z a : we have to show that |n a dj(^) — u(w a )\ > £,a- If 
z G w a w l a for some 1 < i < i(a), then there are three possibilities. First, if w a w' l a C Q \ Q! , then 
u a dj = v £ is affine on the side w a w l a , so the claim is trivial. Second, if w a w l a C Q! and z belongs 
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to the cross Zp associated to the vertex wp = w l a , then again the claim is immediate since u^{z) 
belongs to the ball B(m(^),^), which does not intersect B(u(w a ), £ Q ) by Remark 5.4 Lastly, 
if w a w l a C Q' and z ^ Zp, then u a ^(z) = u(z), thus the claim is again obvious by the definition 
of&. 

To conclude the step, we have to consider a point z ^ Z a which does not belong to any side 
of Q starting at w a . We have again to distinguish some possible cases. If z belongs to the cross 
Zp for some /3, then again the claim follows by the fact that B[u(w a ),^a) n B(u(wp),£,p) = 0. 
If z does not belong to any cross and z G Q' , then Uadj(z) = u(z) so the claim follows because, 
using the bi-Lipschitz property of u and the fact that £ a < £ a /(3L), we have 



u(z) £ B(u(w a ),£ a ) 



\Z — W n \ < 



which is impossible because \z — w a \ > £ a . Finally, consider the case when z £ Q\Q' ■ In this 
case, we surely have \z — w a \ > r by construction, thus by (4.4) and (5.10) we get 

l^adjO*) - u{w a ) \ = \v £ (z) - u(w a ) \ > \u{z) - u(w a )\ - \u(z) - v £ (z)\ 

\Z — W n 



> 



L 



V%r r 
6L3 ~ 2L > ^ ' 



thus the first step is concluded. 



Now, taken two points z, z' £ Q, we have to show the validity of (5.12). 



Step II. Validity of (5.12) if z, z' £ Z a . 

Let us first suppose that both z and z' belong to the same cross Z a . By construction, u a dj 
is L bi-Lipschitz on each segment w a pi, hence to show (5.12) we can assume without loss of 
generality that z £ vo a p\ and z' £ w a p2- Therefore, on one side we have 

l^adjOO ~ u adj(V)| < |«adj(z) ~ «adj(«>a)| + |«adj(™a) ~ ^adj(^')| < L {\ z ~ w »\ + \w a ~ z'\) 



< V2L 



z — z 



On the other side, to estimate |itadj(z) ~~ ^adj^Ol from below, assume without loss of generality 
that liiadj^a) — ^adj(^)| < |«adj(^a) — ^adj (z') | , and define z" £ w a z' so that 

Since the triangle u ai di(w a )u a dj(z)u ax i^(z") is isosceles, then 



7T 



«adj(^adj(z")«adj(^) > 7T ■ 



Moreover, we claim that 



hadj(^) - «adj(^")| 



> 



\Z — Z" 



2L 



(5.13) 



(5.14) 



Indeed, if both w a w\ and w a w^ belong to Q', then by definition 

|"adj(^) - Uadj(z")\ _ |«adj(Pl) ~ Wadjfe) | _ \ u (Pl) ~ U (P2)\ 



> 



z — z" 



Pi - P2 



Pi ~ P2 



L ' 
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so (5.14) holds true. Conversely, if both w Q w^ and w a w 2 a belong to Q\ Q! , then since v £ is L + e 

V £ (P2)\ 



bi-Lipschitz we have 

hadj(^) - 
\z - 



Padj(Pl) 
|P1 



«adj(P2) 



\ve(j>l) 
\Pl 



> 



P2\ 



L + e 



so again (5.14) holds true. Finally, assume that WaW^ C Q' while u? a ii; 2 C Q\Q' (the case 
of w a Wa Q Q\Q! and WaW^ C Q' being completely equivalent). In this case, it must clearly 
be ui Q S <90 e , hence by Remark 5.4 we know that \p\ — w a \ and \p 2 — w a \ are both at least 
£, a /L = r/(3L 2 ), thus \pi — p%\ > \/2r/(3L 2 ). Therefore, recalling again (4.4), we have 



hadj(-z) ~ ^adj(z'' 



> 



|Pl — P2 I 

|u(pi) - u(p 2 ) 

\P1 ~ P2 I 



\u(Pi) - 



«(P2)| 



■P2 



|n(p 2 ) 


- «e(P2)| 




-P2| 



> 



1 \/2r/(6L 3 ) _ 1 
L ~ V2r/(3L 2 ) ~ 2l' 



thus (5.14) has been finally checked in all the possible cases. This inequality, together with (5.13) 



and again with the fact that u ad j is L bi-Lipschitz on the segment z z C w a p2, yields 

2 V 
V2 t 



Padj 



(z) - Uadj (V) I > -^(hadj(^) - W ad j (z") I + I Uadj ~ «twij(^)|) 



> 



2L 



+ 



L 



> A 
^4L |Z 



Summarizing, under the assumptions of this step 

V2, 



4L 



z'| < Uadj(^) 



^adjl 



< y/2L\i 



(5.15) 



Therefore, (5.12) is shown and this step is concluded. 



Step III. Validity of (5.12) if for all a one has z, z' £ mtZ a . 

Consider now the situation when neither z nor z' belong to the interior of any cross. In this 
case, we have that u^(z) = u(z) if z G Q', while u a dj(z) = v £ (z) if z G Q \ Q! , and the same 
holds for z' . Since u is L bi-Lipschitz while v £ is L + e bi-Lipschitz, the validity of (5.12) is 
obvious if both 2, z' £ Q', as well as if both z, z' £ Q \ Q' . Therefore, we can just concentrate 
on the case in which z € Q', z' € Q \ Q'. 

In this case, the main observation is that \z — z' I > V2r/(3L 2 ), since both z and z' must be 
at distance at least r/(3L 2 ) from any vertex w a G dQ, £ , because they do not belong to any cross 
Z a . As a consequence, again by (4.4) we get 

|«adj(z 



«adj( 



> 



\u(z) — v e 
\z — z'\ 



(z')\ > \u(z) - u(z')\ - \u(z) - v £ {z')\ 



V2r \z 
6L 3 " 



2L 



while 



|«adj( 



Madj 



(z')\ = \u(z) — v £ (z')\ < \u(z) — u(z')\ + \u{z) — v £ {z')\ 



< L\z 



V2r 
6U 



< [L + 



2L 



\z — z 
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thus also in this case (5.12) is proven (keep in mind that, since u is a L bi-Lipschitz map, then 

(5.16) 



of course L > l!). In particular, under the assumptions of this step one has 



\z — z | ,3 , 

< \u a dj(z) - UadjOz') < - L \z - z'\ 



2L 



Step IV. Validity of (5.12) if z G Z a and for all f3 one has z' ^ int Zp. 

We pass now to consider the case when z belongs to some cross Z a , while z' does not belong 
to the interior of any cross. In particular, we can assume that z E w a p\. To get the above 
estimate in (5.12), it is enough to make a trivial geometric observation, namely, that there 
exists 1 < i < i(a) such that 

I 'l^^fl l_U 

\z - z\> —\\z ~Pi\ + \Pi - I 



not necessarily with i = 1. As a consequence, we can use the estimate (5.15) of Step II for the 
points z and pi -which both belong to Z a - and the estimate (5.16) of Step III for the points pi 
and z' -none of which belongs to the interior of some Zp— to get 

|«adj(^) - u^(z)\ < |u adj (z) - U ad j(pi)\ + |u adj (pi) - M ad j(z')| < V2L\z-pi\ + - L\p { - z'\ 

< 3 V2L\z-z'\. 
~ 2 1 1 



On the other hand, to get the below estimate in (5.12), let us recall that by Step I we have 

Ua,d)(z) £ Clos B(u(w Q ),(, a ) , '"adj(^) ^ B(u(w a ),^ a ) . (5.17) 

Since u & ^(z) belongs to the radius u{w a )u a ^{p\), then an immediate geometric argument 
from (5.17) implies, as already observed in (5.8), that 

| / \ , A | . \u&d]{z) - U ad j(pi)| + |«adj(Pl) - ^adj(^)| 

KdjOO - u ad j{z )\ > J '-^ ■ 



(5.18) 



Thus, using the L bi-Lipschitz property of tt a dj in the segment w a p\, and the estimate (5.16) of 
Step III for pi and z', we get 



PadjCz) - Uad){Z )\ > — h 



3L 6L 

Summarizing, under the assumptions of this step we have 



z'\ \z-z'\ 

> 



6L 



\z — z 
6L 



< l^adj^) - -U-adj^')! < 7^ V2L\z - z'\ 



(5.19) 



hence in particular (5.12) is again checked. 



Step V. Validity of (5.12) if z G Z a and z' G Zr. 

The last situation which is left to consider is when z and z' belong to two different crosses. 
This situation will be very similar to that of Step IV. Indeed, for the above estimate in (5.12) 
we can again start observing that for some 1 < i < i(a) it must be 
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Then, we use the estimate (5.15) of Step II for the points z, p, t E Z a , and the estimate (5.19) 



of Step IV for the points z' £ Zp and pi -which does not belong to the interior of any cross- 
getting 

- ^adj(V)| < hadj(^) - «adj(Pi)| + |«adj(Pi) ~ "adjOOl 

3 

< V2L\z-pi\ + - V2£|p, - z'\ < 3L\z - z'\ . 



Finally, to find the below estimate in ( |5.12[ ) we notice again that (5.18) holds true, and we use 



the L bi-Lipschitz property of it ac jj in w a pi and the estimate (5.19) of Step IV for p\ and z' , 
obtaining 



|«adj(^) - 'UadjU )| > — h 



z\ \z — z 
> 



3L 18L ~ 18L 

Thus, we have finally checked ( |5.12| ) in all the possible cases and the proof is concluded. □ 

6. Finitely piecewise affine approximation on polygonal domains 



In this last section we give a proof of Theorem |1.2| In fact, the proof is quite short, since it 
is just a simple adaptation of the arguments of Section [5} 



Proof of Theorem 1.2: First of all, assume that SI is an r-right polygon and that u is r-piecewise 
affine on SSI according to Definition 3.5 Then, as already underlined in Remark 5.7 we can 



slightly modify the proofs of Proposition 4.1 and Proposition 5.1 to get what follows. First 
of all, there exist some r such that f G rN, an r-right polygon Sl e CC SI which is part of 
the r-tiling of 0, and an L + e bi-Lipschitz and piecewise affine function v £ : Q, e — > M 2 for 



which (4.1), (4.2), (4.3) and (4.4) hold. Moreover, there exists also a finitely piecewise affine 
map v £ : f2\O e — > A\A e which is C±L 4 bi-Lipschitz and which coincides with u on <9f2 and with 
v £ on d£l £ . Therefore, gluing v £ and v £ exactly as in the proof of Theorem |1.1[ we immediately 
get the required C\L A bi-Lipschitz and (finitely) piecewise affine approximation of u. 

Consider now the general situation of a polygon Q with a map u which is piecewise affine 
on d£l. Of course, there exist a right polygon f2 and a (finitely) piecewise affine and bi-Lipschitz 
map $ : SI — > 0, having bi-Lipschitz constant C = C{Q). The map uo$ _1 is a CL bi-Lipschitz 
map from the right polygon Q to A, which is piecewise affine on the boundary. Then, we can 
apply the first part of the proof to get an approximation v : S7 — > A which is finitely piecewise 
affine and C\C L 4 bi-Lipschitz. Finally, : SI — > A is a C\C 5 L bi-Lipschitz approximation 
of u as desired. Thus, the proof is concluded by setting C'(S7) = C 5 . □ 

Remark 6.1. Observe that the (best) constant C'{0) depends on the geometric features of Q, 
such as the minimum and the maximum angles of its boundary. However, by the construction 
above one has that C'(Q) = 1 whenever Q is a right polygon. 
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